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ABSTRACT 

We discuss the recent discovery by Oppenheimer et al (2001) of old, cool white dwarf 
stars, which may be the first direct detection of Galactic halo dark matter. We argue 
here that the contribution of more mundane white dwarfs of the stellar halo and 
thick disk would contribute sufficiently to explain the new high velocity white dwarfs 
without invoking putative white dwarfs of the dark halo. This by no means rules out 
that the dark matter has been found, but it does constrain the overall contribution 
by white dwarfs brighter than My w 16 to significantly less than 1% of the Galactic 
dark matter. This work confirms a similar study by Reyle et al (2001). 

Key words: Galaxy - dark matter; Galaxy - structure 



> 

CN 
CN 
O 
CI 

o 

6 ' 

in 
c3 



1 GALACTIC DARK MATTER AS FAINT 
STARS? 

The MACHO and EROS microlensing projects have de- 
tected a microlensing signal from dark objects in the Galac- 
tic halo (Alcock et al. (2000) (The MACHO collaboration), 
and Lassarre et al (2000) (The EROS Collaboration). Taken 
together, their results show the microlensing can be ex- 
plained by a Galactic dark halo 20% of which is in the form 
of approximately 0.5 Mq objects. This suggests that the re- 
sponsible objects could be low mass main sequence red dwarf 
stars, or white dwarf stars, but in order to have escaped de- 
tection to date, they must be very faint. 

Red dwarfs are too luminous, or they would have been 
detected directly in the Hubble Deep Field (HDF) (Flynn 
et al 1996, Elson et al 1996 and Mendez et al 1996). White 
dwarfs were considered more difficult to rule out directly, 
until the surprising discovery was made that white dwarfs 
which have had long enough to cool, cease becoming fainter 
and redder, but remain at approximately constant luminos- 
ity while becomeing bluer, due to the development of H2 in 
their atmospheres which induce very non-blackbody spectra 
(Hansen, 1999a, 1999b). 

The possibility that very faint, blue objects had been 
found arose when the HDF was imaged at a second epoch 
and analysed by Ibata et al (2000). The number of moving 
objects, their colours, magnitudes and proper motions were 
all consistent with the detection of a significant fraction of 
the dark halo matter in the form of old, cool white dwarfs. 
However, a third epoch observation of the HDF did not con- 
firm the proper motions of the objects (Richer, 2001), but 
the idea that the dark matter had been detected had already 
spurred many groups to search for local counterparts — with 
some degree of success. As a result, a number of very low lu- 



minosity white dwarfs have turned up in new proper motion 
studies by Ibata et al (2000), Hodgkin et al (2000), Scholz et 
al (2000), Goldman (2000), de Jong et al (2000) and Oppen- 
heimer et al (2001), along with a very low luminosity white 
dwarf identified by Ruiz et al (1995), now viewed with new 
significance. A colour magnitude diagram for these objects 
in shown in Figure [l]. The recently discovered white dwarfs 
are shown as triangles — they are all fainter than the end 
of the white dwarf cooling sequence at Mr « 15.5, and have 
velocities typical of the spheroid white dwarfs. 

Oppenheimer et al (2001 — hereafter OHDHS) have 
conducted the largest of these recent surveys. They discov- 
ered 38 high velocity white dwarfs and derived a space den- 
sity for their objects which corresponds to approximately 2% 
of the Galactic dark halo density at the Sun (aproximately 
0.01 Mq pc -3 ). This is a small but significant fraction of the 
dark halo density, not high enough to explain the microlens- 
ing events (which require that about 20% of the dark halo 
be in the form of ~ 0.5 Mq objects, but still significantly 
higher than the expected contribution of white dwarfs from 
all the well understood Galactic stellar populations — disk, 
thick disk, and stellar halo (or spheroid). 

A very similar survey to OHDHS, in terms of the local 
volume of space surveyed for high proper motion stars, is 
the Luyten Half Second catalog, Luyten (1979). Flynn et al 
(2001) have searched the LHS and two other older proper 
motion surveys for nearby dark halo white dwarf candidates. 
The LHS covers more than half the sky, has a limiting mag- 
nitude of V = 18.4, with proper motions in the range 0.2 to 
0.5 arcsec/year. A recent independent analysis (Monet et al 
2000) shows that the LHS is substantially (90%) complete 
within these limits, based on a new, deeper survey over a 
small area within the LHS. 
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Table 1. Adopted Luminosity function for disk M dwarfs 
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Table 2. Adopted Luminosity function for disk white dwarfs 
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Figure 1. Colour (B — R) versus absolute magnitude (in the 
ij-band) diagram of nearby white dwarfs. Disk white dwarfs are 
shown by crosses, and spheroid white dwarfs as circles. Until re- 
cently, the faintest white dwarfs known were at Mr ~ 16, but re- 
cently four white dwarfs have turned up in new surveys designed 
to probe for them (squares). These objects are shown as triangles. 
One object, J0050-5152, is a binary, the secondary being much 
fainter and its status as a white dwarf is still to be confirmed. All 
of these new white dwarfs are still relatively bright, and probably 
not faint enough to be good dark matter candidates. 



Gibson and Flynn (2001) have searched the LHS for 
objects of the type detected by OHDHS, but found at most 
a few, even though the survey covers a similar volume. If 
2% of the local dark matter is composed of white dwarfs, 
then a few tens of objects had been expected in the LHS, 
whereas Flynn et al (2001) had earlier analysed the LHS and 
two other proper motion surveys in detail, and had found 
no convincing evidence that any of the high proper motion 
objects could be associated with "dark halo white dwarfs". 
All the objects were broadly consistent with coming from 
the visible Galactic populations. 

An alternative explanation for the OHDHS white 
dwarfs is that they are from existing Galactic populations, 
and do not represent a new population from the dark halo. 
That they are members of the thick disk has been argued 
by Reid et al (2001) and Hansen (2001). We also argue for 
this view in this paper by modelling the expected numbers of 
white dwarfs which would be found in proper motion surveys 
from the known Galactic populations. Our work supports a 
similar, independent study by Reyle et al (2001). 



2 MODEL OF LOW LUMINOSITY STARS IN 
THE SOLAR NEIGHBOURHOOD 

We have built a model of the low mass stellar content of 
the Solar neighbourhood, including the luminosity function, 
density and kinematics, in order to simulate actual proper 
motion surveys via a Monte-Carlo technique. The model al- 
lows us to predict the expected number of low mass stars 
which would be recovered in a proper motion survey directed 
toward any point on the sky, covering a given area, with a 
given apparent magnitude limit and a detection window of 
proper motions. For example, the LHS catalog covers a little 
more than half the sky centered on the Northern hemisphere, 
has an apparent magnitude limit of R = 18.6, and recovered 
proper motions (j, in the window 0.5 < (X < 2.5 arcseconds 
per year. 

We model a sphere centered on the Sun with a radius of 
up to 288 pc, which is sufficiently distant from the Sun for all 
the surveys we consider. Only low luminosity (My > 12.5) 
stars are included in the model. 

2.1 Populations in the Model 

The local Galactic components represented in the model are 
the disk, the thick disk, the stellar halo and the dark halo. 

• The disk component consists of M dwarfs and white 
dwarfs. The M dwarfs are drawn from a luminosity function 
shown in Table [l], which has been measured from faint star 
counts with HST (Zheng et al 2001). The white dwarfs are 
drawn from the luminosity function shown in Table ^ which 
comes from Liebert et al (1988). 

The disk white dwarfs have velocity dispersion compo- 
nents a — (au , ay , aw) (where U, V, W are the usual space 
velocities in the directions of the Galactic center, Galac- 
tic rotation and perpendicular to the Galactic plane) of 
a = (40, 30, 20) km s~ 1 , and a mean motion (asymetric drift) 
relative to the Sun of —20 kms . For the disk M dwarfs, 
we adopt velocity dispersion components of a = (35, 25, 15) 
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kms -1 , and a asymetric drift, V aBB = 15 kms -1 . These val- 
ues are slightly lower than for the white dwarfs, because 
the mean age of the M dwarfs is certainly lower than the 
white dwarfs. The young M dwarf component is measured 
to have a = (30, 17, 12) kms -1 , and the old M dwarf compo- 
nent a = (56, 34, 31) kms -1 , in an analysis of the Hipparcos 
results by Upgren et al (1997). In both the surveys anal- 
ysed in this paper, the number of M dwarfs recovered from 
the proper motion windows is quite sensitive to the kine- 
matics adopted; our values have been choosen to produce 
about the right number of M dwarfs. Note though thata our 
main interest is to predict the number of white dwarfs in 
the proper motion surveys, while the M dwarfs are only in- 
cluded as a consistency check on the modelling. The number 
of white dwarfs predicted in the surveys is less sensitive to 
the adopted kinematics, because the white dwarfs are gen- 
erally closer than the M dwarfs. 

• The thick disk component of the model consists of white 
dwarfs selected from the same luminosity function as the 
disk white dwarfs, but with space density of 5% of the disk 
(Reyle and Robin, 2001). This is consistent with recent es- 
timates of the thick disk density as between 2% and 10% 
of the disk density (Kerber et al 2001). We adopt an un- 
certainty in the thick disk local normalisation of about a 
factor of two. The thick disk stars are given velocity disper- 
sion components a = (80,60,40) kms -1 , and an asymetric 
drift V&BB = 40 kms -1 (see e.g. Morrison et al 1990). We also 
include thick disk M dwarfs, selected from the same lumi- 
nosity function as the disk M dwarfs, with a normalisation 
of 5% of the disk, and with the same kinematic parameters 
as the thick disk white dwarfs. 

The significantly larger space motions of the thick disk 
stars means that they are several times more likely to be 
found in a typical proper motion survey than disk stars. 
Reyle et al (2001) model the thick disk in their study of 
white dwarf proper motions as a — (67,51,42) kms -1 and 
Vkss = 53 kms -1 . We have found that adopting either ours 
or Reyle et al's kinematics for the thick disk ends up giving 
quite similar results, i.e. the model predictions are not very 
sensitive to the adopted thick disk kinematics. 

• The stellar halo, (or spheroid) part of the model con- 
sists of white dwarfs drawn from a luminosity function due 
to Liebert (2001, private communication). The luminosity 
function is based on 7 white dwarfs, identified as members 
of the "halo" on the basis of a high tangential velocity, 
Vtan > 160 kms -1 obtained as part of an ongoing analy- 
sis of all the faint stars in the LHS. The luminosity function 
is shown in Table ^| Note that we have converted bolometric 
magnitudes to V-band magnitudes using Eqn 1 of Liebert 
et al (1988). The total number density in these objects is 
3.2 x 10 -5 stars pc -3 which corresponds to a mass density 
of 1.9 x 10 -5 M Q pc -3 for a white dwarf mass of 0.6 M . 
This is approximately 15% of the local mass density of the 
stellar halo (as seen in subdwarfs) of 1.5 x 10 -4 M Q pc -3 
(Fuchs and Jahreifi, 1997). 

For the stellar halo we adopt velocity dispersion compo- 
nents of a = (141, 106,94) kms -1 (Chiba and Beers, 2001) 
and an asymetric drift of 180 kms -1 (i.e. there is a small 
net rotation in the Galactocentric coordinate frame). The 
adopted values have very little impact on the conclusions of 
the paper, because stars with such high velocities populate 
the proper motion selection window quite well. For exam- 
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Table 3. Adopted Luminosity function for stellar halo white 
dwarfs 
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Table 4. Kinematic parameters of populations in the model 
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pie, changing these values to a — (131,106,85) kms -1 and 
an asymetric drift of 229 kms -1 (as used by Reyle et al 2001) 
only changes the predicted number of white dwarfs by a few 
percent. 

We further include stellar halo M dwarfs, drawn from the 
disk M dwarf luminosity function, but with a local density 
reduced by a factor of 500 (see e.g. Morrison 1993). They 
are assigned the same kinematics as the white dwarf stellar 
halo. 

• The dark halo part of the model consists of white 
dwarfs only. We adopt velocity dispersion components of 
a = (156, 156, 156) kms -1 and an asymetric drift of Vas = 
220 kms -1 , which simulates an isothermal population with 
a density falloff with Galactic radius of p oc R~ 2 , i.e. a suffi- 
cient but not necessary condition to explain the flat Galac- 
tic rotation curve. Our starting point for the mass density 
of these stars is 2% of the local dark halo density of 0.008 
Mq pc -3 (Gates et al 1998) with an average WD mass of 
0.6 Mq , as adopted in the OHDHS survey. 

We have adopted a very simple luminosity functions for 
the dark halo white dwarfs. We give all the dark halo white 
dwarfs the same absolute magnitude, Mv = 15.9, which is 
the faintest absolute magnitude of any of the putative dark 
halo white dwarfs found in the OHDHS survey. We show 
that this choice, combined with a 2% dark halo, leads to 
significantly larger numbers of dark halo white dwarfs in 
both surveys than actually observed. Any other choice of 
luminosity function which is still consistent with the WDs 
in the OHDHS sample (i.e. all stars are brighter than My = 
15.9) would produce an even greater disagreement between 
the simulations and the LHS and the number of recovered 
WDs in the OHDHS sample. 

The dark halo white dwarfs share similar kinematic 
properties with the stellar halo white dwarfs (i.e. high ran- 
dom velocities) , but differ from them in a key respect worth 
pointing out (although it has no effect on the modelling or 
conclusion in this paper). The density distribution of the 
stellar halo is well determined by luminous stars, and fol- 
lows a power law outward from the Galactic center, p(R) oc 
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R~ , where R is the Galactocentric radius (see e.g. Wet- 
terer and McGraw 1996). If the dark halo were made en- 
tirely of dim white dwarfs, they would require a density dis- 
tribution which follows p(R) oc R~ 2 , in order to generate 
a gravitational field which would account for the observed 
fiat rotation curves of disk galaxies. However, we know from 
the microlensing results that they do not dominate the to- 
tal mass of the halo. In that case, they could still have an 
R~ S 5 distrubution (as discussed by Gates and Gyuk 2001). 
In their model, "normal" Cold Dark Matter makes up the 
rest of the mass distribution and does follow an R~ 2 dis- 
tribution. For the purposes of the modelling here, the local 
volume surveyed is so small that the density of the white 
dwarfs over the volume is essentially constant. 

The kinematic properties of all the model populations 
are shown in Table H. 



3 MODELLING THE LHS AND OHDHS 
SURVEYS 

The model allows us to run Monte-Carlo simulations of a 
proper motion survey. Stars are generated within a small 
sphere around the Sun, with a uniform density distribution 
in the case of the halo and thick disk populations, and with 
a density distribution which is falling off in z as sech 2 (z/z^) 
in the case of the disk (where we adopt Zh = 125 pc, which 
is equivalent to an exponentially falling disk scale height of 
250 pc far from the plane). 

For each star a V-band absolute magnitude is selected 
from the appropriate luminosity function, and its apparent 
magnitude computed. A V — I colour is assigned to the stars 
as follows: 

For the disk M dwarfs we use the empirically calibrated 
relation due to Reid (1991) 

(V-I)md = 0.297 x My - 0.858 

For the white dwarfs, V — I colors are calculated simi- 
larly by a relation: 

{V - I)wd = 0.385 x M v - 4.85 

This relation comes from fitting white dwarfs in the 
sample of Bergeron et al (2001). We have also derived trans- 
formation equations for white dwarfs from V — I to other 
colours from the same sample: 

B- V = 1.05 x (V- /) - 0.13 

R- 1 = 0.473 x (F- -0 + 0.01 

Note that a small Gaussian random number with stan- 
dard deviation a — 0.05 is added to the V — I colour to 
avoid clutter in the figures. 

An absolute magnitude Mv versus colour V — I diagram 
from a typical simulation of the LHS catalog is shown in 
Figure ^. Main sequence disk stars dominate the sequence 
on the right side. The left sequence consists of white dwarfs 
of the disk, the thick disk and the stellar halo. 

3.1 Simulation of the OHDHS survey 

The OHDHS survey covers 4165 square degrees, centered 
on the South Galactic Pole, to a limiting magnitude of 
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Figure 2. The main components of the model in the Mv versus 
V—I colour magnitude diagram. A small random number is added 
after assigning the color to avoid stars settling strictly on a line. 
The main sequence stars and white dwarfs are clearly separated 
in this plane. 



R59F = 19.7 in a proper motion /1 window of 0.33 < fj, < 3.0 
arcsecs/year. In order to simulate the OHDHS sample, we 
generate colour magnitude diagrams of nearby stars us- 
ing the model described in the previous section. We then 
transform the Mv and V — I values to the bands used in 
the OHDHS survey. For white dwarfs, we use the following 
transformations to the (photographic) R59F and the (pho- 
tographic) Bj bands (i.e. for the Ill-aJ emulsion) used in 
the OHDHS survey: 

R - R59F = 0.006 + 0.059 x (R - I) - 0.112 x (R - I) 2 

-0.0238 x (R-I) 3 

B-Bj = 0.28 x(B-V) for - 0.1 < B - V < 1.6 

which come from Bessell (1985) and Blair and Gilmore 
(1982). 

For the M dwarfs we have derived a transformation from 
R-I to Bj-R 59F based on Figure 10 of Hambly et al (2001), 
who describe in detail the Super COSMOS sky survey, upon 
which the OHDHS survey is based. 

Figure ^ shows the results of a simulation in the reduced 
proper motion Hr versus colour B.j — R plane (hereafter, Hr 
and the _R-band refers to the R59F band, used by OHDHS). 
Panel (a) shows the simulation, and panel (b) shows the ac- 
tual OHDHS data. Circles mark disk stars, squares mark 
thick disk stars and triangles stellar halo stars. There is 
good agreement between the simulated sample and the ac- 
tual data. Firstly, the disk white dwarfs form a wide se- 
quence from (Bj - R, H R ) = (0.0, 17.0) to (Bj - R, H R ) = 
(1.5, 22.0). Most white dwarfs appear to be members of the 
disk. Below this sequence, most of the thick disk and stel- 
lar halo white dwarfs appear, because they have generally 
greater space velocities than the disk stars, and thus higher 
reduced proper motions. 
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Figure 3. Reduced proper motion versus colour plots for the 
OHDHS and simulated samples. Panel (a) shows the data from 
OHDHS, with white dwarfs marked by stars and M dwarfs by 
dots. Panel (b) shows a typical simulation of the OHDHS survey 
using our model. The main features of the data are well repro- 
duced by the simulations: most of the white dwarfs detected are 
plausibly from the disk, thick disk and stellar halo. Two of the 
OHDHS white dwarfs (at high reduced proper motion) are not 
easily accounted for by the model, and are very interesting as 
dark matter candidates. 



The M dwarfs lie in an almost vertical line at Bj — 
R ~ 2.0 in both panels. The total number of M dwarfs in 
the simulation is similar to the observations, as a result of 
adjusting the kinematic parameters of the disk to acheive 
this consistency. The number of disk M dwarfs turns out to 
be quite sensitive to the disk kinematics, and we acheived 
this good fit by using the kinematic parameters for the disk 
shown in Table ^. The kinematic parameters adopted are 
however quite consistent with the observed kinematics of M 
dwarfs (Upgren et al (1997), and also section 2.1). 

To a limiting magnitude of R = 19.7, OHDHS identify 
a total of 97 white dwarfs, of which 82 have direct spec- 
troscopic confirmation, and 15 are assumed with reasonable 
confidence to be white dwarfs based on their position in the 
reduced proper motion versus colour diagram. 

We have computed the expected numbers of various 
white dwarf types using the model. The expected numbers 
are: for the disk 60 ± 8 WDs, for the thick disk 21 ± 5 Wds 
and for the stellar halo, 10 ± 3 WDs, for a total of 91 WDs. 

The adopted kinematical parameters have a small ef- 
fect on the predicted numbers of white dwarfs. For example, 
adopting the disk, thick disk and halo kinematics used by 
Reyle et al (2001), (i.e. for the disk: a = (42.1,27.2,17.2) 
kms -1 , and V as = 16.6 kms -1 , for the thick disk a = 
(67,51,42) kms - , and V as = 53 kms -1 , and for the halo 
a = (131, 106,85) kms -1 , and V as = 229 kms -1 , compare 
with Table §), we obtain 54 disk WDs, 18 thick disk WDs 
and 11 halo WDs, that is the numbers are very similar. 

The total number of white dwarfs in the model and the 
OHDHS sample are in good agreement. We now check the 
relative numbers of white dwarfs of each population type. 
The simulations indicate that a neat dividing line between 
disk and other types of white dwarfs can be drawn from 
(Bj - R, H R ) = (-0.5, 17.0) to {Bj - 7?, H R ) = (2.0, 26.0). 
Counting white dwarfs below this line in the simulations 
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Figure 4. Kinematics of the white dwarfs. Panel (a) shows the V 
versus the U velocities for the white dwarfs in the Oppenheimer et 
al sample. Panel (b) shows results from a typical simulation using 
disk, thick disk and stellar halo white dwarfs, but no dark halo 
white dwarfs. (The same symbols are used as in Figure |jj(b)). The 
main features of the observations appear to be well reproduced 
by the classical stellar populations alone. 



yields 27 ± 5 stars, compared to 21 white dwarfs in the 
OHDHS sample. 

Oppenheimer et al (2001) argue that many or all of 
these higher reduced proper motion white dwarfs are part 
of a new dark halo population, comprising some 2% of the 
local dark matter halo density. We argue instead that these 
white dwarfs can just as well be interpreted as coming from a 
mixture of the thick disk and stellar halo. Further evidence 
for this assertion comes from the kinematics of the white 
dwarfs. 

We show in Figure ^ the space velocities of the white 
dwarfs in the V versus U plane, i.e. projected onto the Galac- 
tic plane. A typical simulation is shown in panel (a) and the 
OHDHS sample stars in panel (b). In the simulation, disk 
stars are marked by diamonds thick disk stars by crosses and 
halo stars by circles. The similarity between the diagrams is 
striking. In particular, the thick disk and stellar halo white 
dwarfs in the simulation lie mostly outside the 2-cr circles 
used by Oppenheimer et al (2001) to isolate the dark halo 
stars (squares in squares). (The 2-cr circle is where stars with 
velocities twice that of the disk velocity dispersion lie, rela- 
tive to the mean motion of the old disk at (V, U) = (—35, 0) 
kms -1 ). In a typical simulation we count 41±6 white dwarfs 
outside the 2-cr circle, while OHDHS find 37 such stars. The 
distribution of these simulated stars in the (V, U) plane is 
also found to be very similar to the observations. 

We show in Figure | theF versus U velocities for a sam- 
ple of local dwarf stars for which [Fe/H] is available (Fuchs 
and Jahreifi, private communication 2001). The availability 
of abundances means that the stars can be classified by their 
population type directly, rather than compared statistically 
to the models (as is the case for the white dwarfs). There are 
282 stars in total, of which 245 are in the metallicity range 
— 1.6 < [Fe/H] < —0.5, and are here termed "thick disk" 
stars, while 37 have [Fe/H] < —1.6 and are here termed 
"halo" stars. What is striking about this figure is how sim- 
ilarly the (V, U) distribution of the stars is to the OHDHS 
samlpe (Figure ^(a)). Thick disk stars dominate the region 
around (V, U) = (—35, 0) kms -1 ), while the halo stars fill a 
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Figure 5. Kinematics of a sample of nearby dwarfs (Fuchs and 
JahreiB, private communication 2002). Crosses indicate stars in 
the range —1.6 < [Fe/H] < —0.5 (broadly, "thick disk" stars) 
while circles are stars with [Fc/H] < —1.6 (broadly, "halo" stars). 
The distribution of velocities for the two metallicity types is very 
similar to that of the white dwarfs of the OHDHS sample (Figure 
|(a). 



much broader region centered roughly at (V, U) = ( — 220, 0) 
kms^ 1 ). 



3.2 Simulation of the LHS 

The LHS covers about 28,000 square degrees, mostly of the 
Northern sky, to a apparent magnitude limit of R = 18.6 and 
with a proper motion window of0.5</i<2.5 arcsec/year. 
We simulate the LHS by covering the fraction of the sky 
(5 > -33° and |6| > 10°, i.e. 65% of the sky) which is 
estimated to be complete to better than 90% by Dawson 
(1986). 

Figure ^| shows the reduced proper motion Hr (RPM) 
versus colour for the simulated LHS sample, where the 
Luyten R band magnitude (Rl) is related to the V mag- 
nitude via (see Flynn et al 2001, appendix B) 

R L = V- 0.17 - 0.228 x (B - V). 

Note that we have not shown the M dwarfs in the simu- 
lation, because we were unable to obtain satisfactory trans- 
formations from V — I to Luyten's photographic Bl — Rl 
colour. 

Although the LHS and the OHDHS surveys probe sim- 
ilar volumes of space, it is much easier to compare our sim- 
ulations with OHDHS. This is because the OHDHS sources 
have been spectroscopically classified into white dwarfs and 
red dwarfs (such work is underway with the LHS and should 
be available in the near future) and furthermore the colour 
transformations are better understood than in the LHS. 
Hence, we regard the comparison of our model with OHDHS 
in the previous section as superior to any comparison with 
the LHS. Nevertheless, the broad features of the LHS ob- 
servations are well reproduced by the simulation, and we 
regard this as a qualitative consistency check on the model. 



Figure 6. Reduced proper motion versus colour for the LHS sam- 
ple. Panel (a) shows the data from the LHS, while panel (b) shows 
a typical simulation (of the white dwarfs only). The white dwarf 
part of the LHS data is broadly reproduced by the simulation. 



3.3 Dark Halo White Dwarfs? 

We now introduce dark halo white dwarfs into the simula- 
tions. We begin by adopting a luminosity function in which 
all the dark halo white dwarfs are as faint as the faintest 
white dwarf found in the OHDHS sample, i.e. at absolute 
magnitude My = 15.9. This choice minimises the number of 
dark halo white dwarfs generated in the simulations, while 
not being inconsistent with the luminosities of the white 
dwarfs actually found in the OHDHS survey. Still, it appears 
to produce too may simulated dark halo sources. Following 
OHDHS, we adopt a white dwarf mass of 0.6 Mq and local 



density of 2% of the dark halo, i.e. 0.02 x 0.008 

-3 



1.6 x 10"" 



Mq pc _il . In Figure (?| we show where these dark halo white 
dwarfs lie in the reduced proper motion versus colour plane. 
Despite being conservative and minimising the number of 
"dark" white dwarfs, there are plenty of these white dwarfs 
at a reduced proper motion of Hr ~ 25, where just a few 
stars are found in the observed samples. 

Judging from the simulations, a good discriminator be- 
tween white dwarfs of the disk, thick disk and stellar halo 
and those of the dark halo is to divide them at Hr = 24. 
In the OHDHS sample, two white dwarfs are found with re- 
duced proper motion Hr > 24, while none were found in the 
LHS. Both samples survey very similiar local volumes. We 
plot in Figure || the number of high reduced proper motion 
white dwarfs (Hr > 24) expected in the LHS and OHDHS 
surveys as a function of their fraction of the dark halo. As- 
suming that no dark halo white dwarf candidates were found 
in the LHS, and two (the two absolutely faintest stars) were 
found in Oppenheimer et al's survey, we conclude that ap- 
proximately 0.4% of the dark matter density could be in 
white dwarfs (i.e. 2 WDs were found corresponding to a 
dark matter fraction of ~ 0.1% in figure]^). A dark matter 
fraction of 2%, as suggested by OHDHS, would yield about 
12 high reduced proper motion white dwarfs in the OHDHS 
sample and about 5 in the LHS sample, for a total of 17, 
compared to 2 WDs actually found. Judging from both sur- 
veys, we conclude that OHDHS 's estimated fraction of 2% 
of the dark halo density in white dwarfs is an overestimate, 
and should be approximately an order of magnitude smaller. 
If a significant fraction (> 0.1%) of the dark halo is in white 
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Figure 7. Reduced proper motion versus colour plots for the 
OHDHS and simulated samples. Panel (a) shows the data from 
OHDHS, with white dwarfs marked by stars. Panel (b) shows 
a typical simulation of the OHDHS survey using our model, now 
including white dwarfs from the dark halo (the figure is otherwise 
the same as Figure []). Two percent of the dark halo mass has 
been placed in white dwarfs of 0.6 Mq at an absolute magnitude 
of My = 15.9, and these are shown in the simulated survey by 
filled triangles. Most of the white dwarfs in OHDHS are plausibly 
from the disk, thick disk and stellar halo, but two high reduced 
proper motion white dwarfs in OHDHS do lie in the region where 
dark halo white dwarfs with an absolute magnitude of My = 
15.9 are expected to lie, although (as shown in Figure ^) the 
simulations produce approximately a factor of ten more such stars 
than actually seen in the combined LHS and OHDHS surveys. 




.005 .01 .015 

Fraction of local Dark Halo Mass 

Figure 8. Expected number of high reduced proper motion 
(Hr > 24) white dwarfs in the two surveys analysed, as a function 
of the fraction of the local dark matter density in such objects 
(for a white dwarf mass of Mwd = 0.6Mq. The actual number 
of objects found in both surveys with this high a reduced proper 
motion is 2, which corresponds to a local density of fa 0.2% of the 
dark halo. The dark halo density adopted is 0.008 Mq pc~ 3 



dwarfs, then they should be fainter than absolute magnitude 
Mv w 16 in order to avoid being detected in large numbers 
in either the LHS or OHDHS. 



We have built a model of the kinematics and luminsities 
of low mass stars in the Solar neighbourhood, which we use 
to simulate the results of various proper motion surveys. We 
discuss in particular the search for very faint white dwarfs in 
the Luyten Half Second proper motion survey and recently 
by Oppenheimer et al (2001), which are the largest surveys 
of the type. The surveys sample similar volumes of space. 
We argue that the contribution of "normal" white dwarfs of 
the thick disk and stellar halo is sufficient to explain the high 
velocity white dwarfs found by OHDHS, and that they are 
not necessarily part of a new, massive population from the 
dark halo. This work confirms results from a similar study 
by Reyle et al (2001). 
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